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Abstract 

. In the classical Monge-Kantorovich problem, the transportation 

cost only depends on the amount of mass sent from sources to des- 
tinations and not on the paths followed by this mass. Thus, it does 
not allow for congestion effects. Using the notion of traffic intensity, 
we propose a variant taking into account congestion. This leads to 
an optimization problem posed on a set of probability measures on a 
suitable paths space. We establish existence of minimizers and give 
s * . a characterization. As an application, we obtain existence and varia- 

tional characterization of equilibria of Wardrop type in a continuous 
■ space setting. 
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1 Introduction 

Given two mass distributions /io and \L\ on IR d with equal total mass, the 
classical Monge-Kantorovich problem consists in finding transportation plans 
(i.e. measures on M. d x IR d having /i and fii as marginals) with minimal aver- 
age transportation cost. This old problem can be traced back to Monge [TTj . 
It has received a lot of attention in the recent years since the pathbreaking 
paper of Brenier [I] who solved the case of a quadratic transportation cost. 
We refer to the book of Villani [12], the lecture notes of Ambrosio pQ and 
the references therein for a recent account of this rich mathematical theory 
and its numerous applications. 

An interesting case is when the transportation cost is given by a confor- 
mally flat Riemannian distance: 



where n(/i ,/ii) is the set of transportation plans. However, from a traffic 
planning point of view, problem ( 11. 1ft is not totally realistic. On the one 
hand, ( II. ip is path-independent: the total transportation cost only depends 
on the amount of mass transported from the sources x to the destinations 
y and not on the paths followed by this mass. Put differently, in (11.11) . in- 
dividual's travelling strategies are irrelevant. On the other hand, (11.11) does 
not take into account congestion effects i.e. the possibility that the cost g(x) 
for passing through the point x depends on "how crowded" small neighbour- 
hoods of x are. This idea can be made precise thanks to the notion of traffic 
intensity associated to a probability measure on a suitable set of paths. This 
notion of traffic intensity (see paragraph [212]) is the path-dependent analogue 
of the well-known notion of transport density in Monge's problem (see Bou- 
chitte, Buttazzo and Seppecher [3], Bouchitte and Buttazzo [2J, Caffarelli, 
Feldman and McCann [6]). With this notion at hand, we propose an optimal 
transportation problem with congestion. This variant of (11.11) takes the form 
of a relatively simple convex optimization problem posed on a set of proba- 
bility measures on a suitable path space. We obtain existence of minimizers 
(theorem I2.9p and a characterization (Theorem 13.81) . 

Researchers in the field of applied traffic modelling have long emphasized 
the role of congestion in networks. In the early 50's, Wardrop (see |13j ) 




the corresponding Monge-Kantorovich problem then reads as 
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considered the situation where a large number of vehicles have to go from 
one location to another, connected by a finite number of different roads. 
Each vehicle has to choose one road (or a probability on the set of possible 
roads) to minimize some transportation cost which depends not only on the 
road chosen but also on the total flow of vehicles on this road. Wardrop 
gave a minimal stability requirement for transportation strategies: the cost 
of every actually used road should be equal or less than that which would be 
experienced by a single vehicle on any unused road. This natural equilibrium 
concept has been very popular since its introduction because of applications 
to networks of course but also due to the development of non-cooperative 
game theory in the meanwhile. From the best of our knowledge, the study of 
Wardrop equilibria have mainly been restricted to the case where admissible 
roads are given by a finite graph. A secondary contribution of the present 
paper is to introduce an equilibrium concept of Wardrop type in a continuous 
state setting, to prove the existence of such equilibria and to relate it to the 
optimal transportation problem with congestion (theorem 14.21) . 

2 Optimal transportation with congestion 
2.1 Notations 

Given a locally compact separable metric space X, we will denote respectively 
by M. + (X) and Ai+(X) the set of positive and finite Radon measures on X 
and the set of Radon probability measures on X. If X and Y are locally 
compact separable metric spaces, /x G A4\(X), and / : X — > Y is a Borel 
map we shall denote by /Jj/x the push forward of /x through / i.e. the element 
of M\{Y) defined by fMB) = fi(f-\B)) for every Borel subset B of Y. 

In the sequel, C d denotes the <i-dimensional Lebesgue measure. If /x and 
v are in M. 1 + (R d ) then ^ denotes the Radon-Nikodym derivative of /x with 
respect to v. We shall write \i << vto express that /x is absolutely continuous 
with respect to u, in which case, slightly abusing notations, we will identify 
/x with the Radon-Nikodym derivative j^. 

The data of our problem are Q (its closure Q modelling the city, say) 
which is some open bounded convex subset of M 2 , two probability measures, 
/x and H\ in A4\(Q), giving respectively the distribution of residents and 
services in the city Q. The set of transportation plans associated to /xo and 
Hi will be denoted II (/j , /xi) it consists of the probability measures on Q x Q 
having /x and /ii as marginals: 

IT(/xo, /xi) := {7 G M\(Ti x O) : 7T h = )"o, = M (2-1) 
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where (vr (x, y), iri(x, y)) := (x, y), stand for the canonical projections {x and 



Introducing congestion naturally leads to consider spaces of paths, lengths 
of such paths and sets of probability measures on sets of paths. From now, 
on we shall denote: 

• C := W 1 ' 00 ^, 1], 0), viewed as a subset of C°([0, 1],R 2 ), 

• C x ' y :={a EC : a(0) = x, a(l) = y} (x, y in O), 

• 1(a) := J Q \a(t)\ dt, the length of a E C, 

• for a E C, a denotes the arclength reparameterization of a belonging 
to C, hence \cr(t) \ = 1(a) = 1(a) for a.e. t E [0, 1], 

• C := {a E C : \a\ is constant} = [a, a E C}, 

• for Q E M\(C), we define Q E M\(C) as the push forward of Q 
through the map a i— > a, 

• for (p E C°(n, R) and a E C, we define 



• e (a) := a(0), e x (a) := a(l), for all a E C°([0, 1], M 2 ). 
2.2 Traffic congestion modelling 

The classical Monge-Kantorovich optimal transportation problem for a given 
cost function c E C°(Q x f2,R) is: 



Note that, in the linear problem (j2.2p . the cost of transporting one unit of 
mass from x to y, c(x, y), is given and does not depend on the path(s) followed 
by the mass from x to y. In order to take into account congestion effects, 
we explicitely introduce probabilities over C x,y as part of the optimization 
problem. More precisely, the overall transportation cost will depend not 
only on the transportation plan 7 E il(yU ,/ii) but also on the way travelers 
commuting from x to y use the different possible paths a E C X)V . In the 
sequel, the way commuters from x to y are split according to the different 



y in Q). 





(2.2) 
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paths will be given by a probability measure p x ' y on C x,y . Put differently, 
jf' y (E) is the proportion of travelers from x to y using a path a G £ C C^. 
This naturally leads to the following definition: 

Definition 2.1. ^4 transportation strategy consists of a pair (j,p) with 7 G 
Il(/i ,/i 1 ) and where p = ip x ' y )t Xty ) e Uxn ^ s a B or el family of probability mea- 
sures on C such that p x ' y (C x ' y ) = 1 for^-a.e. (x,y) Gflxfi. 

There results, from the use of a transportation strategy (7,p), an overall 
traffic intensity J 7iP G .M+(f2) defined by 

[<p(x)dln(x) ■=[__([ ( f <p(*(tW(t)\dt) dp^(a)) d 7 (x,y) 

= LS\ L Mdp x,y (°)) d i^y)i Vy? G c°(n,K) 

(2.3) 

and an overall probability over paths Q 7iP G M\(C) given by <5 7 , p = p x,y ® / y, 



i.e.: 



F(<7)dQ 7 » =[__([ F(a)dp x > y (o-)) d^y) VF G C°(C,R). 
Jnxn \Jc x >y J 

(2.4) 

One could consider the probability Q 7 , P as if it represented the total number 
of travelers that use a path a G £ given the global transportation strategy 

(t,p)- 

Let us remark that if we set Q := Q ltP G ,M+(C) then J 7iP only depends 
on Q, and can be written as 1 1V = iq G M.+{VL) where %q is defined for every 
Q< M\[C) by: 



/_ 



(p(x)di Q (x)= [ L v {a)dQ{a), V</? G C°(f2, R). (2.5) 



Let us also remark that since L^a) = L ip (a) one has Iq = ig, for all Q G 
A1+(C). Finally, let us note that the total mass of %q is the average length 
with respect to Q: 

i Q (U) = [ l(<j)dQ(<j). (2.6) 
Jc 

if the probability Q is concentrated on injective curves, one could also express 
the measure %q through H 1 — integrals. In this same case, namely 
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Remark 2.2. If a curve a is injective, then one could also write 1(a) = 
7i 1 (cr([0, 1])) and L v (a) = J^^LpdH 1 . Moreover, if for 7— a.e. (x, y) the 
probability p x,y is concentrated on the set of injectives curves from x to y, 
one could also define the measure J 7iP by replacing the integral with respect 
to \a(t)\dt in (12.31) with an integral in dTC 1 . Notice moreover that for every 
Borel subset 4c0 one would have: 

I y>p (A) =[__([ H l {Ana)dp x ' y {a) \ dj{x,y) = [ H l (A D a)dQ y>p (a). 
Jnxn \Jc x -y J Jc 

If we imagine that for each a £ C x,y , the mass of travelers commuting on 
cr is uniformly distributed on a, this means that I ljP (A) represents the cu- 
mulative traffic through the region A. The same formula stays true, under 
no injectivity assumption, if we replace 7i 1 (A fl cr) with Li A (a) and in this 
case the cumulative traffic takes into account the number of times a path a 
passes through the points of A. 

In the sequel, it will be convenient to formulate our optimization problem 
in terms of Q = Q llP rather than in the transportation strategy ( r y,p). To 
that end, we shall use the following: 

Lemma 2.3. Let us define 

Q(fi ,fii) := {<3 7iP : (j,p) transportation strategy} 
then one has 

Q(/x ,A*i) = {QG M\{C) : e^Q = fi , e^Q = /ij. 

Proof. If (jiP) is a transportation strategy then e t!Q 7 ,p = 7T tj7 = Ho, and 
eiH<5 7 , P = v"i(j7 = /ii. Now let Q £ M.\{C) be such that e$$Q = /io, ei$Q = 
Hi. If we define 7 := (eo,ei)Jj<5, we have 7 £ n( / u 0? / i i)- It then follows from 
the disintegration theorem (see [8]) that there exists p = (p x ' y )( Xty ) € Tixn a 
Borel family of probability measures on C such that p x >y(C x,y ) = 1 for 7-a.e. 
(x, y) £ Q x Q and Q = p x,y ® 7. Hence Q = Q 7 , p for a transportation 
strategy (j,p). □ 

At this point, a natural way to model traffic congestion is, for a given 
transportation strategy (j,p), to consider that the transportation cost per 
unit of mass between x and y is given by 

<hA*,v)= [ L Ghp (a)dp x ' y (a) (2.7) 
Jc x -y 
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where Gi is a nonnegative function which depends (in a way that will be 
specified later on) on the traffic intensity J 7iP . The optimal transportation 
with traffic congestion then takes the form (to be compared with the usual 
Monge-Kantorovich problem ( 12. 2ft ): 

inf < / c JyP (x,y)d'y(x,y) : (7, p) transportation strategy > . (2.8) 

linxn J 

Setting Q = Q 1P and using formally (12.51) . we see that the total transporta- 
tion cost in (12.81) can be rewritten as: 



/_ _c^ p (x,y)dj(x,y) = / L G (a)dQ(a) — _G iQ (x)di Q (x). 
Jttxtt Jc Jn 

Hence using lemma 12.31 we can reformulate (12.81) in terms of Q only: 

(2.9) 



inf iJ_G iQ {x) di Q (x) : Q e Q(po, (ii) 



Note that in the definition ( 12.71) . it is required that Gi p is continuous (or at 
least l.s.c.) whereas the form ( 12.91) allows for more general forms of congestion 
through i ^ Gi. From now on, we assume that G has the following local 
form: 

Gi(*) = 9 (Jg(z)\ , (2-10) 

where ^2 is the radon-Nicodym derivative of i with respect to the Lebesgue 
measure and g is a nondecreasing function R + — > M + such that the function 
H defined by H(z) = zg(z) for all z G IR+ is convex and superlinear (i.e. 
lim 2 _ +00 g(z) = +00). 

The optimization problem we shall study now then reads as: 

inf H(i Q ) where H(i) = { ^ H ^ dx if 1 « C ^ (2.11) 
Q£Q(mo,mi) I +°° otherwise. 

In the sequel, we shall say that a transportation strategy (7,2?) is optimal if 
Q 7iP solves (12.131) . 

Remark 2.4. It will be clear in the sequel that the probability Q 7)P associated 
to an optimal transportation strategy (7,^) will be concentrated on injective 
curves, so that the interpretation in terms of Ti 1 — integrals (see Remark 12.2ft 
may apply. 
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2.3 Existence of minimizers 



From now on, we make the following assumptions: 

• H is convex and nondecreasing on R + with H(0) = 0, 

• there exists q > 1, and positive constants a and b such that 
az q < H{z) < b{z q + 1) for all z G M+, 

• H is differentiable on IR + , and there exists a positive constant c such 
that < H'{z) < c{z q ~ l + 1), for all z G R + , 

• the following set 

Q q (fi ^i) ■= {Q e QOo,m) : *q e L"} (2.12) 

is nonempty. 

These assumptions enable us to simply rewrite (12.111) as: 

inf f H(i Q (x))dx. (2.13) 

Remark 2.5. Let us discuss the assumption that Q IJ (/io,A t i) 7^ which, at 
first glance, may seem difficult to check. In order to have the existence of 
a Q G Q(fio,fi\) such that %q G L g it is sufficient that /io and [i\ are in L 9 . 
This result, which is not obvious, follows from the regularity results of De 
Pascale and Pratelli (see [H] and [TO]) who proved that L q regularity of /x 
and fii implies that for 7 solving the Monge-Kantorovich problem (12.21) with 
c(x,y) = \x — y\ and p x,y = 5[ Xty ] (the Dirac mass at the segment [x,y]) for 
every x and y the corresponding traffic density J 7)P is L q . It is not necessary 
however that //o and /ii are absolutely continuous for the assumption to be 
satisfied: let us consider for instance the case where Q = [0, l] 2 and /i and 
fix are respectively the one-dimensional Hausdorf measures of the segments 
[(0,0), (0,1)] and [(1, 0), (1, 1)]. If we define 7 := (id, id + (1,0))^ and 
px,y _ £ then a straightforward computation shows that J 7iP is uniform 
on[0,l] 2 . 

Under the assumptions above, we are going prove that (12.131) admits a 
solution. The proof of existence involves some preliminary lemmas. 

Lemma 2.6. For any if G C (O, M+) ; L v is lower semi- continuous on C for 
the uniform topology, indeed for any a G C , one has: 



H, inf Up o a) J \a(t i+1 ) - a{U)\ : 
1. . \[ti,t i+ i\ 1 

([ti,t i+ i])i is a subdivision of [0, 1] j. (2.14) 
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Proof. For any subdivision ([ij, t i+ i])i = \ y .. n , we have: 

n rU+i 

L M = E / v>M*))l<H*)l dt 
i=l J ti 

n rti+x 

> V inf ( V oa) / |or(t)| dt 
[ti,t i+ i\ j t . 

n 

. \ [U,U+i\ 
1=1 

Taking the supremum over all such divisions, we get: 

n 

L<p(cr) > sup < inf ((poa)\a(ti+i)—a(ti)\ : ([U,t i+1 ])i is a subdivision of [0, 1] 

^ . ' [U,U+i] 



i=l 

Let us prove the converse inequality. Let e > 0, since ip o a is uniformly 
continuous, there is a 5 > such that: 

Vt,t'e [0,1] 2 , (|t-t'| <<5^ -^(0)1 <£)• 

For any subdivision ([tj, £i+i])i=i,... n such that \t{ — U+i\ < 5 for all z, we have: 



fH+l 

L»<]T( inf (<po*)+e) / 

i=1 W.**+lJ Jt, 



dt 



n 



1=1 J 

([Tj,Tj + i])j is a subdivision of 
<sup{VV( inf ((p o <r) + e)\a(Tj) - <r(T j+ i)\ : 

([rj, Tj+i])j is a subdivision of 

n 

= sup{V( inf (pa) + e)\a(t i+1 ) - a{U)\ : 
i=i 

([ti,t i+ i\)i is a subdivision of [0, 1] 



As this last inequality is true for any e > we get (12.141) . The lower semi- 
continuity is then obvious since, by (12. 141) . is the supremum of family of 
lower semi-continuous functions on C°([0, 1], Q). □ 
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Lemma 2.7. Let (Q n ) n G A^(C°([0, 1], M 2 )) N sucA ^ Q n (C) = 1 /or a// 
n and there exists a constant M > snc/i that: 



sup / Z(tr) rfg n (a) < M. 

n Jc 

Then the sequence (Q n ) n is tight and admits a subsequence that converges 
weakly * to a probability Q such that Q(C) = 1. 

Proof. The tightness of (Q n ) n easily follows from the inequality: 
Q n ({aeC: \&\ > #}) = Q n ({a G C : 1(a) > K}) 

< 1 / rf Q n ( a ). (2.15) 

By Prokhorov theorem, we may therefore assume, passing to a subsequence 
if necesseary, that (Q n ) n converges weakly * to Q G M.\(C°([0, 1], M 2 ). It 
remains to show that Q{C) = 1. For > let us define Cx := {cr G 
C : < X}, then Inequality (12. 15ft and the fact that the measures Q n are 
concentrated on C yield 

~ ~ M 
sup Q n (C\C K ) = sup Qn(C\C K ) < — , 

n n A 

for every K > 0, which implies 

1 = limsup<5n(C) < lim sup Q n (C K ) + lim sup Qn(C \ Csr) 

n n n 

M 

< Q(C K )+—. 

Letting K tend to oo, we then get Q(C) = sup K Q(Ck) — 1. □ 

Lemma 2.8. Let (Q n )n be a sequence in Ai + (C) that converges weakly * to 
some Q G M.\{C). If there exists i G A4 + (Q) such that iq n converges weakly 
* to i in A4 + (Q) then we have ig < i. 

Proof. Let y? G C°(U,R + ), we first have: 



/ ipdi = lim / ifdi Qn = lim / L^dQn 
Jn n Jn " n Jc 



it easily follows from lemma [231 that Q \— > j c L^dQ is l.s.c. for the weak * 
topology of J\A\(C), we then have: 

ipdi > / L^dQ = I ipdiQ. □ 
n Jc Jn 
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Now, we are in position to prove: 

Theorem 2.9. The minimization problem A2.13\) admits a solution. 

Proof. Our assumptions imply that the value of (12.131) is finite. Let (Q n ) n 
be some minimizing sequence of (I2.13p . From the identity iq = iq, we may 



assume Q n = Q n for all n. We deduce from our growth condition on H, 
that {iQ n ) n is bounded in L q . On the one hand, extracting a subsequence 
if necessary, we may therefore assume that {iQ n ) n converges weakly in L q to 
some i. On the other hand, since iQ n is bounded in L q and hence in L 1 we 
have 



sup / l(a)dQ n (a) = sup / i Qn 
n Jc n Jn 



< +00. 



Moreover Q n = Q n and we deduce from lemma 12.71 that (up to some sub- 
sequence) (Q n )n weakly * converges to some Q in M. l + {C). Since Q(/io, /ii) 
is obviously weakly * closed, we have Q G Q(/Uo,/ii) and lemma |2~%1 implies 
that %q < i (consequently to this enaquality iq is absolutely continuous). 
From the monotonicity and convexity of H we then have: 

H(iQ(x))dx < / H(i(x))dx < liminf / H(iQ n (x))dx 
Jn n Jn 

which proves that Q solves (12.131) . □ 

Let us remark that if H is furthermore assumed to be strictly convex then 
if Qi and Q 2 solves (12.131) then i Ql = %q 2 so that the optimal traffic intensity 
is unique (of course, this does not imply in general that Qi = Q 2 or that the 
corresponding optimal transportation strategy is unique). 

3 Characterization of the minimizers 

In the sequel, we shall denote by q* the conjugate exponent of q, given by 
q* = q/(q-l). 

3.1 Optimality conditions 

The variational inequalities characterizing solutions of the convex problem 
(I2.13P can be expressed as follows 

Proposition 3.1. Q G Q 9 (Ai ,/i 1 ) solves \2.13\) if and only if 

[ e% = inf{ / Bq ■■ Qe QVc/ii)) withl:=H'(iQ)eL«*. (3.1) 
Jn Un ) 
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Proof. Assume that Q solves (12.131) . then for every Q G <2 g (/i , A^i), one has: 

< lira -\H{iq^q_q)) ~ = Hm -[H{iq + e(* Q - %)) - 

= / #'(%)(*Q-%) = / £0'q-%) 

which proves fl3.ll) . Conversely, if Q £ Q q (no,^i) satisfies fl3.ll) . then by 
convexity of if, for every every Q G Q q (fio, Hi), one has: 

w(»'q) - w(ny) > / - %) > o. □ 

The next paragraphs will be devoted to investigate the precise meaning 
of (13. ip . Before going further, let us do some formal manipulations to give a 
formal interpretation of ( 13. ip in terms of optimal transportation strategy. Let 
us assume that Q solves (I2.13P and let us write Q = Q-^ for some (optimal) 
transportation strategy (7, p) and define £ := H'^ig), then (13. ip formally can 
be rewritten as: 

/ = f L-^a)dQ(a) 
Jn Jc 

L^a)dp^ y (^d-f(x,y) 

inf f (f L^a)^{a))Mx,y) 

isp. strategy JnxTi \Jc*>v J 



nxn \Jc x >y 



( 7 ,p) transp. strategy JnxTi \Jc*< 

inf / inf / Lj{a)dp(o-) d , y(x,y) 

inf / ( inf Lj(a) ) ^7(2, y) 



7Sn(Ai ,A 

defining (again formally) the transportation cost: 

c i (x,y) = Jnf s L ? (a), 

we then firstly have: 

/ c?(x,y)dj(x,y) < / LjdQ = inf / c?(x,y)dj(x,y) 
JnxTi Jc 7en(uo,Mi) ./nxn 

so that 7 solves the Monge-Kantorovich problem: 



inf / cdx,y)d~{(x,y). 
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Secondly: 



Lj(a)dQ(a) = _c^(x,y)d>y(x,y) 

C JSJxSl 



c^(a(0),a(l))dQ(a) 
c 



and since Lj(a) > cj(a(0), a(l)), we get 

Lj(a) = qr(a(0), cr(l)) for Q-a.e. a. 

or, in an equivalent way, for 7-a.e. (x, y) one has: 

L|-(ct) = c^(x,y) for p x ' y -a..e. a. 

Since £ is only L 9 * , and are not well-defined and the previous arguments 
are purely formal. In paragraph !3.2l we will extend the definition of eg to the 
case where £ is only L q * under the additional assumption q < 2. This will 
enable us to make the formal argument above rigorous and to characterize 
optimal transportation strategies in paragraph 13.31 



3.2 The transportation cost when £ is L q * 

For a non-negative function £ G C°(Q) we define 

q c (x, 3 /) = inf{L € ( < 7) : a G C*>»}. 

Proposition 3.2. Let us assume that q < 2 and define a := 1 — 2/q* , then 
there exists a non-negative constant C such that for every £ G C°(fi,M + ) and 
every (xi, yi, x 2 , y 2 ) £ one ^ as -' 

| C f(x!, j/0 - c € (x 2 , y 2 )| < C||e|U (n) - x 2 \ a + \ yi - y 2 \ a ) . (3.2) 

Consequently, if (£ n ) n G C°(f2,R + ) N bounded in L q * , then (c^ n ) n admits a 
subsequence that converges in C°(Q x f2,R + ). 

Proof. Let £ G C°(n,R+) and ar, y G fi 2 . For jfe > let a k G C x > y be such 
that 

1 1 

t{vk{t))\a k {t)\ dt<cs(x,y) + -. 

Then for all e > and z such that y + ez G and to G (0, 1) we consider the 
following element of C x ' y+ez : 



t7 fc (^) ifte[0,t ] 
f(^)ez iftG[t ,l]. 
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We then have, for all k > 0: 

ct=(x, y + ez) < I i(a kM {t))\a kM {t)\ dt 



Z(a k (t))\a k (t)\ dt + / £(y + tez)e\z\dt 
Jo 

f 1 1 

< W, y) +e\z\ J £(y + tez)dt + -. 



Now we let k tend to oo and we get 

c 6 (x,y + ez) -cz(x,y) 

e 

and by a similar argument 



< \z\ [ i{y + tez) dt, 
Jo 



c,(x, y) - c,(x, y + ez)_ ^ f ^ 



e 

r l,oo 



This implies that c%(x, .) G W ,0 ° and: 

\V y ct(x,.)\ < |£(.)|, for alia;. (3.3) 

By symmetry we also have 

\V x ct(.,y)\<\£(.)\, for ally. (3.4) 

Since q* > 2, we deduce from (13.31) . ( 13. 4p and Morrey's Theorem (see [5], 
Chapter IX), that there is a constant C such that: 

\ct(x,yi) - C((x,y 2 )\ < C\\£,\\ Lq * \y x - y 2 \ a , for all x,y h y 2 in fi, 
|c e (xi,?/) - c^(x 2 ,y)\ < C||^|| L ,*|a;i - x 2 \ a , for all x u x 2 ,y in fi. 

This proves (13.21) . The second claim in the proposition then follows from 
(HOD, the identity Cg n (x,x) — and Ascoli's theorem. □ 

From now on, we further assume that q < 2. For a non- negative function 
£ G L 9 *(f2) we then define 

C((x, y) = sup {c(x, y) : c G .4(0} , 

where 

.4(0 = {limc en in C\H x O) : (£„)„ G £ n > 0, £ n -> £ m L 9 *} . 
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Remark 3.3. The definition of is unchanged if one replaces £ ra — > £ in L 9 
by — ^ £ i n £ 9 * i n the definition of A(£). Indeed, if we do so, we obviously 
obtain a function which is larger than eg. Now, let us assume that £ n — i - 
£ in L g * , and c^ n converges to c in C°(f2 x Q), using Mazur's Lemma there 
exists a sequence 7/ n which converges strongly to £ and such that each r^ n is 
in the convex hull of k < n}. It is clear that for fixed x, y, £ — > c^(x,y) 
is concave hence c(x,y) = lim c^ n (x, y) < limsup c Vn (x,y) < c^(x,y). 

When £ is continuous, one has: 

Lemma 3.4. If £ is continuous and non-negative, then — c^. 

Proof. The inequality > is obvious, as one can always choose the con- 
stant sequence £ n = £ in the definition of eg. Let us show now the opposite in- 
equality. Take x, y G f2, e > and a G C 3 ^ such that L^(a) < c^(x, y) + 1/k. 
We can choose cr so that it is piecewise linear, by density of this kind of curves 
and using the continuity of £. Let (<Si)i=o,...,m-i be the segments which com- 
pose a with Si = XiXi + \, xq = x and x rn = y. Let us fix, moreover, a sequence 
£n £ such that — * c. Now, we want to prove c < c^. Fix a small number 
5 > and for any a G [0, 5] let us define a curve a a in the following way: let 
R be the clockwise 90 degrees rotation in the plane; let x'^a) and x"(a) be 
the only points such that Xix'^a) = aRei and Xi + \x"(a) = aRei, where is 
the tangent unit vector to a in the Si part; define a a by linking any point 
x'^a) to x"(a) by some segments S'^a) and x"(a) to by some arcs 

A i+1 (a) with center x i+1 and radius a. In this way we have a a G C Xa ' Va , 
where x a = x' (a) and y a = x'^(a). Let Ri(S) be the rectangle whose vertices 
are the points x'l{8) and and let Bi(8) be the circular sector 

centered at Xi and whose vertices are and rr -(5). 

If we compute L^ n (a a ) da it is not difficult to see that we get 

„S m-l , „ \ m-1 / „ x 

/ L u (a a ) da = J2( indC 2 + £ / ■ 
Moreover it holds c^ n (x a ,y a ) < L^ n (a a ), hence we get 

„ 5 m-l ✓ „ v m-l ✓ „ x 

/ c^ n (x a , y a )da<J2( UdC? + E / ^ ■ 

If we pass to the limit as n — > oo we get, by using the uniform convergence 
of c^ n to c on the left hand side and the L q * convergence of £ n to £ on the 
right hand side, 

/ c(^,y«) da<Y,( Zd£ 2 + £ / ^ 2 . 
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Then we divide by 5 and we pass to the limit as 5 — > 0. Using the fact that 
c is continuous we have 



1 f 5 

V m „ 7 / c ( x °' y a ) da = c ( x ' y)- 



On the other side, we may notice that the areas of the sectors Bi may be 
estimated by CS 2 and hence we have, for 5^0, 



i=i Vfii / 



On the contrary the integrals over Ri, when divided by 5, converge on the 
integrals on the segments Si, which give exactly the integral over the curve 
a, i.e. L^(a). We have consequently 



m—l 




c(x,y) < <%(x,y) + ^ 

and, being arbitrary, we also get c < and the thesis. □ 

Lemma 3.5. Let us assume that q < 2 and let £ 6e non-negative function 
belonging to L q * , then there exists a sequence (£„)„ G C°(Q), £ n > 0, £ n — > 
£ m snc/i that c^ n converges uniformly to c^ on Q. 

Proof. It is easy to see that for every (x, y) G VL 2 there exists a sequence 
of non-negative continuous functions (£ n ) n converging to £ in L 9 *(f2) such 
that c^ n converges in C° and c^(x,y) = lim n c^ n (x, y). Let / be a finite 
set, (xi,yi) G f2 2 for alH G / and for every i, let be a sequence of non- 
negative continuous functions converging to £ in L 9 * (fi) such that c^(xi, yi) = 
lim n c^i (xi, yi). Let us set £ n := max !€ ;^, we then have (£ n )n converging to 
£ in l/* (ft), and 

ctfayi) < liminf c^fo,^) < limsup c in (x i: yi) < c^(x i: yi). 

n n 

We thus have c^(x h yi) = lim n c^ n (xj, y») for every i G J. Now, let (xi,^)i eN 
be a dense sequence of points of Q 2 . From what preceeds, for every n, there 
exists a continuous non-negative such that 

llfn-fllw* < -> 1^(^,2/*;) - c in (x k ,y k )\ < -, VA; < n. 

16 



By the Holder estimate of proposition 13.21 and Ascoli's theorem, passing to 
a subsequence if necessary, we may assume that c^ n converges in C° to some 
c. Since obviously c(xk,yt) = c^(xk,yk) for all k , we deduce c = and the 
desired result follows. □ 

The next lemma enables us to extend in some sense when £ > is 
only L q : 

Lemma 3.6. Let us assume that q < 2. Let Q G Q q (fiQ, fi\) , £ be a non- 
negative element of L g * , and (£ n ) n be a sequence of non-negative continuous 
functions that converges to £ in L q * , then we have the following: 

(i) (£f n ) n converges strongly in L X (C,Q) to some limit which is indepen- 
dent of the approximating sequence (£ n ) n and which will again be de- 
noted Z/£. 

(ii) The following equality holds: 

[ £{x)i Q {x) dx = [ L^a) dQ(a). (3.5) 
Jq Jc 

(Hi) The following inequality holds for Q-a.e. a G C : 

L^(a)>c^(a(0),a(l)). (3.6) 

Proof. For all n and m in N we have: 



c Jc 



(Zn(<r(t))-U(<r(t)))\cr(t)\ dt 



o 



dQ(a) 



< / \£, n (x) - i m {.x)\i Q {x) dx 
Jn 



This implies that (L^ n ) n is a Cauchy sequence in L l (C, Q) and it is obvious, 
from the previous inequality, that its L X (C, Q) limit does not depend on the 
approximating sequence (£„)„. 

The proof of (ii) follows from (i): 

£(x)iQ(x) dx = lim / £ n (x)ig(x) dx 
a n Jn 



lim / L u (a) dQ(a) 
n Jc 

Lt(a) dQ(a). 



c 
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To prove (iii) we choose an approximating sequence as in lemma 

13.51 and pass to the limit in 

^>)>%>(0),cr(l)). □ 

3.3 Characterization of optimal transport strategies 

In this paragraph, our aim is to make the formal arguments of paragraph 
13.11 rigorous in order to characterize optimal transport strategies. This can 
be done under the additional assumption that H is strictly convex. First, 
we relate the optimality condition (13.1 p to the Monge-Kantorovich problem 
with cost c^: 

Proposition 3.7. Let us assume that q < 2 and that H is strictly convex. 
If Q solves 112.1 Sty and £ := H'(iq) then we have: 



t i G = „^¥ ^£*<? = i nf J ^(x,y)d^(x,y). (3.7) 
Proof. Let us recall that from proposition 13.11 we have: 



f% = inf / £i Q . (3.1 



n 



Let £ be a non-negative element of L q * and let Q G Q q (fi ,fii), using 
Lemma l3~o1 and the definition of Q q (no ) Hi) yields: 

ii Q = I L(.dQ > ! c^(a(0)),a(l))dQ(a) 
Jc Jc 



> inf / ct(x,y)d>y(x,y). 
7en(/i ,Mi) Jnxn 

We then have, for all (6i ? *,(>0: 

n * nf ^ / &Q - i nf J c i (x,y)d^(x,y) (3.9) 
Qee«(Mo,m) Jn 7en(/i ,Mi) Jnxn 

and by a similar argument, for all £ G C°(fi,R + ): 

inf £di Q > inf / c^(x,y)dj(x,y). (3.10) 
QeQ(Mo,Mi) Jn 7en(Mo,Mi) Jnxn 

Let £ G C°(0, M+), and e > 0, for every x and y in f2, there exists 
ap y G C x ' y such that (x,y) h- > cr? ,2/ is measurable (see for instance [7j) and 
by Lemma [3.41 

L 4 ((T^) < c 6 (x, y)+e = c^x, y) + e. (3.11) 
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Let 7 G n(/i ,/ii) and let us define the element of Q(p ,pi): Q e := f^p <8>7, 
we then have: 

= /_ _^«' y ) rf 7(^,Z/) < /_ _c^{x,y)d^{x,y) +e. 
Jn Jnxn Jnxn 

Since 7 and e are abitrary, using (13.10j) we obtain 

inf / 'tdi Q = inf / c $ {x,y)d 1 (x,y), V£ E C°(Tl,R + ). (3.12) 

In what follows, for every p G ,M+(0), we extend p by outside f2. Let 
(Pn)n be a standard mollifying sequence. For n G N* let us consider the 
regularized problem: 



inf / H(p n *t Q ) (3.13) 

The existence of a solution Q n of (13.131) can be obtained by similar arguments 
as in theorem 12.91 (using lemma 12.71 and the fact that the L l norm of p n * iq 
equals the total mass of iq). Proceeding as in proposition 13.11 and defining 
jn ■= Pn*i-Q n , £n ■= H'(j n ), T] n := p„ we have: 

/ H'(p n -kiQ ){p n *iQ ) =/ injn= I Vndiq = i?f I T] n di Q 

(3.14) 

With (JSI2D, we then get: 



Vndiq = / injn = inf c Vn (x,y)d>y(x,y). (3.15) 

Jm? 7en(Mo,Mi) Jnxn 

By convexity of if, we also have: 

/ H(j n )< [ H(p n *iq)< [ p n *H(iQ) (3.16) 
Jr 2 Jr 2 ' Jr. 2 

which implies that j n is bounded in L q . Passing to subsequences, we may 
therefore assume: 

j n - j in L\ i n - £ in £ 9 \ tfc - f in L 9 *. (3.17) 

Since the total mass of iq is the same of j n and j n is bounded in L q (and 
hence in L 1 ), we get a bound on iq (0) and, from lemma l2~771 we may also 
assume: 

Q n A Q in .M+(C), ^ i in _M + (H). (3.18) 
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It is obvious that i — j and lemma |2~51 implies j >ig. With (I3.16P and the 
monotonicity of H, we then get: 

/ H(i Q ) < [ H(j) < liminf / H(j n ) < [ (3.19) 
Jn Jn n Jr 2 Jn 

With the strict convexity of H and the optimality of Q, this also yields 

% = i Q = j G L q and liminf / H( 3n ) = [ (3.20) 



Up to some subsequence (as £ = H'{iq) e L 9 *), this also implies 
H(j n ) - H(iq) - 1(j n - ig) — >• a.e. and in L 1 

and using the strict convexity of H, we deduce that j n converges a.e. to %. 
This implies that £ n converges a.e. to £ = H'{i-q) and that £ = £. With 
Fatou's Lemma and (13.151) . we therefore obtain: 

Jn Jk 2 

< lim inf / £ n j n 



liminf inf / c Vn (x,y)dj(x,y). 
n 7 en( W) ,/Ji) JnxTi 



Using r) n — v £ and remark 13.31 from the uniform convergence of c Vn to a cost 
c < cj, we get 



inf / cj(x,y)d-f(x : y) 
n 7en( Mo ^i) JqxQ 

together with (13. 8p and (13.91) . this completes the proof. □ 

The characterization of optimal transport strategies then reads as: 

Theorem 3.8. Let us assume that q < 2 and that H is strictly convex. A 
transportation strategy (^f,p) is optimal if and only if, setting Q := Q^^ and 
£ := H'{i-Q), one has: 

1. 7 solves the Monge-Kantorovich problem: 

inf / c?(x,y)dj(x,y), (3.21) 
7er% ,Mi) Jnxn 
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2. for Q-a.e. a G C , one has: 

Lj{a)= c^(0),a(l)). (3.22) 

Proof. Let us assume first that the transportation strategy (7, p) is optimal 
and set Q := Q-^ and £ := if'(ig). From Proposition 13.71 and Lemma l3~ol 
we get: 

c^(x,y)dj(x,y) = / ^(<r(0),<r(l))dQ(«7) 
<>><<> ic 

< / L ? rfQ = / 

inf / CF(x,y)dj(x,y) 



7en(/i ,/ii) 



SlxSJ 



this proves that 7 solves (13.211) and implies that the inequalities above are 
equalities. We therefore deduce (13.221) from the inequality C|-(cx(0), a(l)) < 
Lfc) . 

Conversely, assume that the transportation strategy (7,^) satisfies the 
two conditions of the theorem. Condition (13. 22f) firstly yields: 



/ L l d Q= I c i{°'^)^{l))dQ(a) = _c^(x,y)dj(x,y) 
! Jc Jc Jnxn 

Secondly, if Q = Q 7:P G Q q (fio, fii), one has: 

fli Q = [ LjdQ > [ c^a(0),a(l))dQ(a)= [_ 
Jn Jc Jc Jn 

and since 7 solves (I3.2ip , we finally have 

Jn Jn 



Ci(x,y)dj{x,y) 



which, with proposition 13. 11 proves that (j,p) is optimal. 

□ 

Remark 3.9. Let us remark that £ = H'{iq) (with Q solving (12. 13ft ) solves 
the following (dual of (12.131) ) problem: 



sup W(0 - / H*(£(x))dx (3.23) 
?ei/?*,£>o Jn 
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where H* is the Fenchel transform of H and: 



W(&) := inf / cz(x,y)dy(x,y). 
Indeed, it follows from proposition 13.71 and £ = H'(iq) that 
W(£)- / #*(£(*))& = / -#*(£)) = / 

Jf2 Jf7 Jf2 

If £ G L 9 * is non-negative, we deduce from (13.91) and Young's inequality: 

W(0< (&q< [ H(ig)+ [ #*(£) 
Jn Jn Jn 

which proves that £ solves (13.231) . 

Let us see, through an easy exemple, an application of Theorem 13.81 

Exemple 3.10. Suppose that Q contains the two segments A = {0} x [0, 1] and 
B = {1} x [0, 1] and the square S = [0, 1] x [0, 1] which is their convex hull. 
Set Hi = 7i}\—A and [i2 = 7i}\—B and denote by T : A — + C the map that 
associates to every point (0,x) G A the curve T(x) given by T(x)(t) = (t,x), 
i.e. the horizontal segment from A to B starting from x. Set Q = TJj/ii. It 
is clear that Q comes from an admissible transportation strategy linking \i\ 
to fi2 and it is not difficult to see that the traffic intensity %q has constant 
density 1 on S and elsewhere. We consider two particular cases only: we 
claim that Q is optimal if Q =]0, l[x]0, 1[ while it is not if S is compactly 
contained in Q. Indeed, if Q — S, the metric induced by %q is the euclidean 
metric, the paths T(x) are geodesic and the transport plan induced by Q is 
optimal according to this metric. On the other hand, if Q is larger than S, 
then all the segments T(x) that are very close to the upper or lower boundary 
of S are not geodesic accoding to this metric, because they could be improved 
by non-straight line paths which arrive up to zone Q \ S where iq = and 
the trasportation is cheaper. In the former case, consequently, the sufficient 
optimality conditions are satisfied, while in the latter the geodesic conditions 
on the paths (Wardrop condition, see the nect section) is not and prevents 
optimality. 



4 Application to equilibria of Wardrop type 

In this final section, we relate the results of the previous sections to some 
concepts of equilibria of Wardrop type. Modelling congestion as in paragraph 
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12.21 enables us to extend the concept of Wardrop equilibrium to a continuous 
setting. 

Let us consider a congestion function g : R + — > IR + which is continuous 
increasing and satisfies az 9_1 < g < b(z q ^ 1 + 1) for all 2 G M+ and some 
g G (1, 2) and non-negative constants a and 6. Then for any transportation 
strategy (j,p) such that I 7iP (defined by (12.31) ) belongs to L q , the transporta- 
tion cost function resulting from the strategy (j,p) is for £ := goI lp G L 9 *. 
Roughly speaking, an equilibrium is then a transportation strategy (7, p) that 
satisfies Wardrop stability condition (i.e. Q ljP gives full mass to the set of 
"geodesies" for the metric £ = g o J 7p ) and the additional requirement that 
7 is an optimal transportation plan between /x and /ii for the cost resulting 
from (7,2?). This leads to the following 

Definition 4.1. ^4 transportation strategy (7, p) zs said to fre an equilibrium 
if i^p G and, setting £ := g o J-^ one aas 

1. L|-(cr) = C|-(cr(0), a(l)) for Q-^-a.e. a E C, 

2. 7 solves the Monge-Kantorovich problem: 



Only the first condition above is linked to Wardrop's original equilibrium 
concept. Imagine that some social planner chooses the transportation plan 
7, then the second equilibrium condition expresses that 7 is optimal for the 
transportation cost resulting from 7 itself and the traveler's individual be- 
havior. Our notion of equilibrium can therefore be viewed as a refinement 
of the Wardrop equilibrium or its generalization to the case where the trans- 
portation plan is not given a priori. 



A direct application of theorems 12.91 and 13.81 then gives the existence of 
equilibria together with a variational characterization: 

Theorem 4.2. Under the assumptions of this paragraph, there exists an 



Remark 4.3. A slightly different situation, which can be relevant in some 
applications, occurs when the transportation plan 7 G n(/i , fix) is fixed and 
not only the marginals /Jo and fi±. In this case, one defines equilibria as the 




equilibrium. Moreover (7,2?) is an equilibrium if and only if Q := Qj t p solves 
the minimization problem: 




(4.1) 



23 



set of p's such that (7,2?) satisfies the first condition (Wardrop) of definition 
SH If the set: 



Q?(j) : = {Q E M\{C) : (e , e^Q = %i Q E L q } 

is nonempty, then slightly adapting our arguments, we have existence of equi- 
libria and p is an equilibrium if and only if Q := Q^-p solves the minimization 
problem: 
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